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ABSTRACT 

The dynamic response of a droplet vaporization process excited hy 
traveling transverse acoustic oscillations is derived "by linear analysis. 
Results of a previous nonlinear numerical study used to formulate the 
analytical model are reviewed, and a transfer function representative of 
the dynamics of the vaporization process is derived. The analysis provides 
dimensionless parameters related to propellant physical properties that 
characterize the dynamic behavior of the vaporization process. Applica- 
tion is made to the vaporization of heptane, oxygen, fluirine, ammonia, 
and hydrazine. The dynamic response of these propellants attains a peak 
value at a particular frequency. A comparison is made with a "burning- rate 
process described by a characteristic time: and an interaction index giving 
similar behavior. 
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INTRODUCTION 


Acoustic mode instability in a rocket engine combustor occurs when 
the combustion energy is released in a manner that reinforces the acoustic 
oscillations. Various combustion processes that can limit or control 
energy release (chemical kinetics, drop burning, vaporization, jet breakup, 
drop shattering, etc. ) have been suggested as responsible for acoustic 
reinforcement. Many of these processes have been studied in a variety 
of investigations. For example, the dynamic behavior of several indi- 
vidual processes is reported in Refs. 1, 2, and 3, and process behavior 
measured in terms of combustor stability is given in Refs. 4, 5, and 6. 

Such studies have isolated specific problem areas in rocket engine insta- 
bility requiring additional study. 

Propellant vaporization is a process of particular interest. Ref. 4 
shows that variations in vaporization rate can significantly affect 
combustor stability. The vaporization process alone was studied more 
thoroughly in Ref. 7 in which the dynamic behavior of the process was 
obtained from a nonlinear numerical analysis by using a model that predicts 
the vaporization rate during the entire drop lifetime. Under certain 
conditions, the process exhibited dynamic behavior that may cause com- 
bustor instability. Such dynamic behavior was related to boundary con- 
ditions controlling the vaporization process of heptane drops. 
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The nonlinear numerical analysis of Ref. 4 motivated the development 
of a linear dynamic analysis "based on a simplified model for vaporization 
that gives the same dynamic "behavior observed in Ref. 7. The advantage 
of the linear analysis is that it provides dimensionless parameters rela- 
ted to propellant properties that may be used to characterize and examine 
the dynamic behavior of the vaporization process for any propellant. 

In this paper the results of the nonlinear numerical analysis of 
Ref. 7 are reviewed, and the linear analysis is presented and discussed. 
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NCMENCLATURE 

nozzle throat area, in. 2 

acoustic velocity at nozzle throat, in. /s 

vapor pressure - liquid temperature coefficient, P^/T^ 

constants 

specific heat of liquid, Btu/(lb)(°R) 
molecular diffusion coefficient, in. 2/s 
frequency, cps 

gravitational constant, 32.2 (lb mass) (ft) /(lb force) (s 2 ) 

thermal conductivity, Btu/ ( in. ) ( s) ( °R) 

mass of propellant being vaporized, lb 

molecular weight, lb mass/(lb) (mole) 

response factor 
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Nusselt number for heat transfer 
Nusselt number for mass transfer 
interaction index 
pressure, lb/in. 2 

combustion chamber pressure, lb/in. 2 

vapor pressure at propellant surface, lb /in. 2 

Prandtl number 

heat transfer rate, arbitrary units 

heat transfer to propellant, Btu/s 

heat transfer from propellant surface, Btu/s 

universal gas constant, 18 510 in. -lb/(°F) (lb) (mole) 

initial drop radius, p 

Reynolds number 

instantaneous drop radius, in. 

Schmidt number 

Laplace transform, d/dt 

temperature of vapor film, °R 

combustion gas temperature, °R 

temperature of propellant, °R 

static temperature at nozzle throat, °R 

time, s 
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period of oscillation, l/f, s 
t 5 Q time to vaporize 50$ of drop mass, s 

final combustion gas velocity, ft/s 
volume, arbitrary units 
vaporization rate, lb/s 
nozzle mass flow rate, lb/s 
correction factor for heat transfer 
correction factor for mass transfer 

vapor pressure - combustion chamber pressure parameter, 
dimensionless 
ratio of specific heats 
gas density at nozzle throat, lb/in. ^ 
latent heat of vaporization, Btu/lb 
characteristic time, s 
mean drop lifetime, M/w, s 
phase shift, deg 
frequency, rad/s 

All primed quantities denote perturbation quantities (i.e. , x' = (x - x)/x ) 
and all barred quantities denote mean values 
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RESPONSE FACTOR 


A response factor can tie defined that is one measure of the magni- 
tude "by which the combustion process can reinforce an acoustic oscilla- 
tion. Such a response factor was introduced in the numerical analysis 
of Ref. 7 and will he used in this study to evaluate the dynamic behavior 
of the vaporization process. This response factor is based on the 
Rayleigh criterion for acoustic amplification by heat or mass addition. 
The Rayleigh criterion states that reinforcement or amplification occurs 
when an excess of heat or mass is added while the pressure is greater 
than the mean value. The response factor is defined as the integral 
value of such excess heat or mass addition over a given period of time 
in a finite volume normalized by the magnitude of the pressure perturba- 
tion. It is an effective gain or amplification factor and is expressed 
for perturbations about a mean value by the following relation 


N = 


rr q’(t,V)P’(t,V) 
f />(*,*■ <a\ 


dV dt 


( 1 ) 


where q' is the fractional heat or mass perturbation, 
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and P' is the fractional pressure perturbation 


P’ 



P 


When "both q' and P' oscillate with the same periodicity and are 
uniform over a finite volume, the response factor is usually expressed 
for one period of oscillation as follows: 



For sinusoidal oscillations in pressure, any heat release or mass 
flow process linearly related to pressure gives the following value of 
the integral: 



P« 

-unax 


cos 0 


(3) 


where 

P'(t) = P^ax sin wt 

and 

q T (t) = q^ ax Sin (wt + 0) 

In analyses of propellant vaporization, the heat release rate q' 
is generally assumed synonymous with the mass release rate w’ . This was 
assumed in the nonlinear analysis of Ref. 7. Also, an average value for 
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the fractional perturbation in mass-release rate over a finite volume 
was determined "by numerical integration"^. For these evaluations, w* 
varies nonlinearly with sinusoidal oscillations in pressure. The response 
factor as defined by Eq. (2) , however, was approximated by numerical 
techniques over one period of the pressure oscillation by the integral 



This integral gave nearly identical values of N to that of Eq. (2) 
for the periodic oscillations in w' observed in Ref. 7. Eqs. (2) and (4) 
give exactly identical values of N only for linear behavior. 

Numerical values of the response factor within the range of -1 to 1 
were obtained in Ref. 7. Negative values indicate that excess heat is 
added when the pressure is less than the mean pressure (potential damping 
of acoustic oscillations), and positive values indicate that heat is 
added at pressures above the mean pressure (potential driving of acoustic 
oscillations) . 

Some significance can be placed on actual numerical values of the 
response factor if the acoustic system is assumed to consist only of a 
heat or mass addition from vaporization (a potential acoustic gain) and 
a mass loss through an exhaust nozzle (a potential acoustic loss). If 
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quasi- steady behavior is assumed, the combustor pressure and flow pertur- 
bations in a critical flow nozzle are in phase. The magnitude of these 
perturbations for adiabatic flow can be derived from 

where a t = -/ri RT^ giving 



If mass flow perturbations in the nozzle are considered analogous to heat 
or mass flow perturbations in the chamber, then for sinusoidal oscillations 
in pressure, the response factor as given by Eq. (3) is 

H= ( '°' 9l2) ^ 1 - 2 

In a simple feedback loop corresponding to the assumed two-process 
system, the sum of such a negative nozzle response factor and the response 
factor of the vaporization process indicates whether an excess of mass is 
added when the pressure is higher or lower than the mean pressure, and 
thus, whether acoustic oscillations will decay or grow according to the 
Rayleigh criterion. 

The numerical result from this analysis of nozzle flow process 
suggests that the response factor for the vaporization process must 
exceed 0.912 to cause combustion instability. The exact dynamic behavior 
of nozzle flow and the magnitude of other acoustic losses and gains in 
rocket combustors is debatable and, therefore, a precise criterion for 
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stability cannot be established. In this study, therefore, attention 
will he confined to how the vaporization response factor varies with 
propellant properties and combustor operating conditions. 

NUMERICAL ANALYSIS 

The response factor as a function of oscillatory frequency was 
evaluated in the numerical analysis of Ref. 7 for a representative vapor- 
ization process and a specific acoustic mode of resonance. Calculations 
were made for n-heptane drops vaporizing in a cylindrical combustor con- 
taining heptane- oxygen combustion products. Pressure, velocity, and 
temperature oscillations associated with the first traveling transverse 
acoustic mode were superimposed on the normal combustion flow process. 
Drops of constant size were repetitively injected from positions uniformly 
distributed across the injector face. The drop vaporization theory 
developed in Ref. 8 was used. In these calculations, the acoustic oscil- 
lations affected the heat and mass transfer processes and the drop accel- 
eration, which gave three-dimensional velocity components to the drops 
and caused perturbations in drop temperature and vaporization rate. 

The perturbation in vaporization rate at a given angular and radial 
position in which the pressure was uniform at an instant of time along 
the axis of the chamber was obtained from a summation of the vaporization 
rates of individual drops. A typical plot of the vaporization rate per- 
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turbation for times covering one complete pressure oscillation is shown 
in Fig. 1. Eq. (4) was used to evaluate the response factor for each 
angular and radial position. The response factor for the entire chamber 
was taken as the average value of the response factor for all angular 
and radial positions. 

These calculations were made for a range of oscillatory frequencies 
and for a variety of boundary conditions affecting drop vaporization 
including variations in combustor pressure, final combustion gas velocity, 
drop radius, initial drop temperature and velocity, and the amplitude of 
the pressure perturbation. The frequency response curves calculated for 
these boundary conditions were correlated by a frequency parameter as 
shown in Fig. 2. An equally good single curve representation of results 
was obtained from a correlation based on drop lifetime. This curve is 
given in Fig. 3, which shows response factor as a function of a dimension- 
less time equal to the ratio of the drop half-lifetime (time to vaporize 
one-half the droplet mass) to the period of the oscillation. The response 
factors shown in Figs. 2 and 3 are larger than those reported?. The 
values have been recalculated and represent the average values of the 
response factors based on the local pressure in each volume rather than 
the wall pressure used in Ref. 7. 


11 



Figures 2 and 3 show the characteristic response of the vaporization 
process for heptane. For fixed boundary conditions, the response factor 
is a •maximum at a specific frequency (corresponding to a dimensionless 
time of about 0. 4) . The response factor approaches zero at a lower 
frequency and a constant negative value at higher frequencies. The tran- 
sition to negative response occurs at a dimensionless time of about 1.4. 
From these response characteristics, the vaporization conditions conducive 
to driving and damping of acoustic oscillations can be readily established. 
It should be noted, however, that the maximum response factor (about 0.8) 
is less than the nozzle loss (about 0.912) previously calculated for 
simple dynamic behavior of exhaust flow. 

The factors responsible for this behavior of the vaporization process 
are of fundamental interest. A description of the contributing factors 
was given in Ref. 7. Parameter groupings that can characterize specific 
behavior, however, are not readily deduced from numerical analysis. Such 
groupings could be readily evaluated if linear dynamic analysis were 
applicable to the problem being studied. 

The primary restriction on linear analysis is the dependency of the 
vaporization process on the absolute velocity difference between the drop 
and the combustion gases. Acoustic particle velocity is a maximum at both 
high and low pressures during the oscillation. For this reason, velocity 
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difference and vaporization rate attain maximum values at conditions near 
loth maximum and minimum pressure and cause the nonlinear "behavior of 
vaporization rate shown in Fig. 1. An analysis of the perturbation curves 
in vaporization rate showed that the velocity difference contribution to 
the vaporization rate at high and low pressure were nearly equal and thus 
cancelled effects with regard to response factor evaluations. As an 
approximation, therefore, the velocity difference effects on vaporization 
rate can be ignored. With this assumption, a linearization of the vapor- 
ization equations is possible. The following section presents a linear 
dynamic analysis of a simplified vaporization process that is insensitive 
to perturbations in gas velocity. 

LINEAR ANALYSIS 

A transfer function that characterizes the dynamic behavior of the 
vaporization process of various propellants will be derived by linear 
analysis of a simplified vaporization model. 

Vaporization Model 

Analytical relations for a vaporizing drop presented in Ref. 8 will 
be used to develop linear equations for perturbations in the vaporization 
process of a single drop. The equations are assumed to apply to a drop 
in an array of drops of decreasing size down a combustion chamber. 
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Drop mass: 


dM 

dt 


= -w 


- dM’ 

Tv dF 


= -w 


where 


M 


t v = • 


Vaporization rate: 


where 


and 




_ (mean drop lifetime) 
w 


2«D/)fr 

w = — — HumPLa 

RT 


= 2 + 0.6 (Sc) 1 / 3 (Re) 1 / 2 


Pc , Pc 

a = wr- 7n 


*5 ‘ U pT^£ 


(5a) 

(5b) 


(6a) 


(6b) 


(6c) 


For this analysis, is assumed proportional to (rP c ) p /^, D is 

proportional to P£ p , and M is proportional to r 3 , so that 

w = CpM in — — 


P c " P L 


and 


w' = 1/2 M’ + pP£ - (p - l/2)P£ 


( 6e) 
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where 


Vapor pressure: 


where 



In P L = C 2 - 


c 3 

t l - c 4 


P L = bT L 


(7a) 

(7h) 


h = 


(T L - C 4 )' 


Drop temperature: 


dT- 


L 1 


dt CpM 


('Jin ~ 'lout) 


(8a) 


If the variations in specific heat Cp with drop temperature per- 


turbations are assumed negligible, then 


CpTL- dTf, 

T V q in ■ q out 


(8b) 


where 


q in - q out “ 
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Heat transfer to drop: 


q in = 2nkNu h r(T b - T L )Z (9a) 

where 

Nu h =2+0.6 (Pr) 1 / 3 (Re) 1 / 2 (gh) 

If the temperature difference (T g - T L ) and the correction factor 

for simultaneous heat and mass transfer Z are constant, and Nu n is 

l/2 

proportional to (rP c ) ' , and M is proportional to r 3 , then 

q in = CsM 1 / 2 !^/ 2 (9c) 

and 

q! n = 1/2 M’ + 1/2 P* (9d) 

Heat transfer from drop : 

lout = Aw ( 10a ) 

If variations in the heat of vaporization A with drop temperature 
perturbations are assumed negligible, 

q’ ut = w' (10b) 

The model expressed by the linear equations may be described as 
follows. A drop with a mass M and a mean lifetime t v is being 
vaporized. The mass of the drop is perturbating about a steady- state 
value because the vaporization rate, which represents mass leaving the 
drop, is varying about a steady-state level. The variations in vapori- 
zation rate depend on the mass of the drop and on the difference between 
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the vapor pressure at the drop surface and the combustor pressure. The 
drop temperature determines the vapor pressure with the temperature 
established in a heat reservoir. The heat to the reservoir depends on 
the combustor pressure and on the mass of the drop or mass in the reservoir. 
Heat leaving the drop or reservoir is contained in the vaporizing propellant. 


Transfer Function 

The linear Eqs. ((5b), (6e), (7b), (8b), (9d) and (10b)) can be 
combined by the use of the transform s = d/dt to give the following 
transfer function involving the variables w' and Pi 


t - Vs [ 

c I 1 + 2 t v s 


1 + (1 - 2|3) ~ 

A(3b 

1 + 2£±t T V s 
Apb 


TyS 


( 11 ) 


This transfer function expresses the dynamic relation between w' 
and P£. Solutions of the form w'(t) = w^ ax sin (cut + 6) can be obtained 
when P'(t) = P-^ax s * n The amplitude ratio and phase shift obtained 

from such solutions are 
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and 


0 = - tan" 1 2 tvW + tan -1 (l - 2|3) — ' £ — Tyto - tan" 1 la. 1 * — Ty<*> 


Apt 


Apb 


( 13 ) 


The response factor defined hy these solutions is given by the 
previously defined relation, Eq. (3), 


” = (!) 

' c 'max 


cos 0 


Propellant Evaluations 

An evaluation of the response factor N for a heptane drop is shown 
in Fig. 4 as a function of the dimensionless time parameter t v co. 

Response factors for oxygen, fluorine, ammonia, and hydrazine together 
with that for heptane are given in Fig. 5. Physical properties used for 
the evaluations are given in table I. The listed properties are for the 
equilibrium drop temperature condition attained during steady vaporization 
at a combustor pressure of 300 psi as given in Ref. 8. For comparison of 
the effects of physical properties on response, Fig. 6 shows the nonlinear 
evaluation of N for the one condition of oxygen vaporization considered 
in Ref. 7. 
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DISCUSSION 


The response curve of heptane vaporization obtained by linear anal- 
ysis (Fig. 4) has characteristic properties similar to that obtained by 
numerical analysis (Fig. 3). Differences in the response curves are 
primarily confined to the region of peak response. The linear analysis 
assumes that vaporization is initiated at an equilibrium drop temperature. 
A numerical evaluation in Ref. 7 for drops introduced at their equilibrium 
temperature gave a peak value of N equal to 0.42, which is the value 
obtained by linear analysis. The 800° R initial drop temperature cal- 
culation shown in Fig. 3 also gives comparable values. Heating of the 
drop from an initial injection temperature apparently introduces nonlin- 
earity that increases the peak value of the response factor. 

A comparison of time bases for the response curves obtained by 
linear analysis of a single drop and by nonlinear analysis of a complete 
array of drops down the chamber shows that t 5 Q for this array is larger 
by a factor of about 4.5 than t v for the single drop. Since t$o is 
the half- lifetime of the largest or initial drop in the array, this means 
that the effective mean drop size of an array of drops is significantly 
less than the largest or initial drop. The relation between the largest 
drop and the effective mean drop depends on the distribution of drop 
sizes in the array. When Eqs. (5a) and (6d) are used to specify time 
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histories, however, an effective mean drop size of about one-fourth the 
initial drop size is obtained. Beyond this limitation of specifying an 
effective drop size, the linear analysis adequately describes the dimen- 
sionless times at which peak, zero, and negative response are obtained. 
These times are characterized by the parameter groupings appearing in 
the linearly derived transfer function, Eq. (ll). Fig. 7 displays the 
contribution to the overall response of heptane vaporization to these 
parameter groupings. 

The term 2x v s/(.l + 2x v s) characterizes the dynamics related to the 
quantity of propellant being vaporized. It establishes the response at 
low frequencies or small vaporization times. In this region, the quantity 
of propellant in the chamber varies inversely with the vaporization rate. 
In the extreme condition of zero dimensionless time, the propellant 
vaporizes as rapidly as it is introduced, and the response factor is 
zero. At frequencies equal to l/2x v and greater, the variations in the 
quantity of propellant being vaporized become small and do not restrict 
the vaporization rate. 

The parameter CpTuv/APb characterizes the dynamics of heat storage 
within the propellant. A frequency equal to l/(cpTLX v AP' b ) essentially 
separates regions of oscillating and constant drop temperature. At lower 
frequencies the drop temperature oscillations approach an equilibrium 
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condition with the oscillating environment. At higher frequencies, the 
temperature oscillations are reduced and eventually drop temperature and 
vapor pressure remain constant. 

A condition of constant vapor pressure gives the inverse effect of 
pressure oscillations on vaporization rate from that obtained when vapor 
pressure is in equilibrium with the environment. Eq. (6e) expresses 
this characteristic. With equilibrium conditions (generally low frequen- 
cies) , an increase in combustor pressure increases the heat transfer to 
the drop, and the vapor pressure and the vaporization rate give a positive 
response factor. When vapor pressure is constant (generally high fre- 
quencies) , an increase in combustor pressure suppresses vaporization 
rate, and a decrease in pressure accelerates the rate, a. condition com- 
parable to flash vaporization. This process gives a negative response. 

The parameter (l - 23) ( CpTLT v /^f5b) characterizes the dynamics of these 
two interacting effects on vaporization rate. The values of 3 for the 
propellants evaluated in this study are greater than unity. Under such 
conditions, the transition from a vapor pressure to a combustion pressure 
controlled process occurs at a frequency that is lower by the factor 
l /(23 - 1) than the frequency characterizing the drop temperature behavior. 

These parameters, which characterize the dynamics of the vaporization 
process, are a function of propellant physical properties. From a knowl- 
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edge of such properties, the dynamic response of various propellants can 
he estimated or computed. Oxygen vaporization (Fig. 5) is one example. 

The region of positive response is broader, and the response factor 
attains a larger value than for heptane. The numerical evaluation (Fig. 6), 
although incomplete, implies a similar behavior. 

Oxygen properties differ from those of heptane in several ways, 
however, the change in (3 (1.36 for heptane; 4.143 for oxygen) is the 

predominant factor causing a difference in response. The value of (3 
for oxygen is large because the equilibrium vapor pressure for oxygen is 
a larger fraction of the combustor pressure than for heptane. When |3 
is large, the dimensionless time characterizing a constant vapor pressure 
increases. This phenomena increases and broadens the region of positive 
response. A large |3 also increases the perturbation in vaporization 
rate when vapor pressure is constant and thus increases the magnitudes 
of negative response as shown in Fig. 5. 

The comparison of various propellants (Fig. 5) shows the region of 
positive response for hydrazine to extend over a broader region of dimen- 
sionless time than for heptane or oxygen. Fluorine is nearly identical 
to oxygen. This comparison implies that a combustion process controlled 
by hydrazine vaporization would be unstable for a broad range of combus- 
tor designs. Ref. 8 indicates, however, that decomposition of hydrazine 
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could effect the vaporization process, which would also effect its 
dynamic "behavior. 

The dynamic response of other propellants for which physical 
properties are known or can he estimated can he evaluated in a similar 
manner. The linear analysis provides a convenient method of surveying 
the dynamic response of propellants and indicates the conditions in 
which the more precise numerical evaluation would he useful. 

STABILITY CRITERIA 

The dynamic response of the vaporization process does not, in itself, 
specify combustor stability. Stability criteria are acquired from a 
complete dynamic analysis of the overall combustion system. Such criteria 
have been derived in Ref. 9 for any combustion process that can be 
characterized by the flowing burning rate expression 

w' = n(p’(t) - P£(t - t)J (14=) 

In this expression, the interaction index n and the characteristic 
time x specify a particular combustion process. The dynamic response 
for this burning rate expression is 

N = n(l - cos x«) (15) 

A comparison of this dynamic response with that of heptane vaporization 
is shown in Fig. 8. The region of positive response and particularly 
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the peak response are of primary significance in establishing system 
stability limits. In this region, the dynamic response of the vaporiza- 
tion process can be approximated by the burning rate expression (Eq. (14)). 
For the comparison shown in Fig. 8, the vaporization process is charac- 
terized by n = l/2 Nmax and t = 4. 5 t v . By a similar comparison with 
the correlation of numerical results shown in Figs. 2 and 3, the value 
of t is approximated by the relation 



within the range of boundary conditions considered in the numerical 
analysis. The value of n is again equal to one-half the peak value 
of the response factor, n ^ 0. 4. If the dynamic response of the vapori- 
zation process is characterized in this manner, the system stability 
criteria specified in terms of n and t are directly applicable to a 
vaporization controlled burning rate. For greater precision, a burning 
rate expression may be postulated that more closely approximates the 
vaporization process and utilizes the method of analysis of Ref. 9 to 
establish stability criterion. 


CONCLUDING REMARKS 

The results of this study show that the dynamic response of the 
vaporization process can be approximated by a linear analysis that assumes 
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the vaporization rate to he insensitive to gas velocity perturbations. 

The assumption is based on the net effect of velocity rather than actual 
insensitivity to velocity, i.e., the changes in vaporization rate due to 
perturbation in velocity are the same at high and low pressures and, 
therefore, cancel with regard to dynamic response. The assumption 
appears valid for any transverse mode in the absence of any steady angu- 
lar or radial flow within the cavity. With such steady flow, however, 
and in the case of longitudinal modes with axial flow of combustion gases, 
dynamic response is not insensitive to velocity perturbations. 10 The linear 
analysis and the results of the numerical analysis do not necessarily 
apply to such conditions. 
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TABLE I. - PROPELLANT PHYSICAL PROPERTIES AT EQUILIBRIUM DROP 
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Fig. 1. - Nonlinear response of vaporization 
process (ref. 7). 



Fig. 2. - Correlation of response factor with vaporization parameters for heptane from 
Ref. 7. Response factors are based on average rather than maximum pressure 
amplitude reported in Ref. 7. 
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Fig. 3. - Correlation of response factor with dimensionless time for heptane from Ref. 7. 
Response factors based on average rather than maximum pressure amplitude reported 
in Ref. 7. 



Fig. 4. - Linear dynamic response of heptane vaporization. 
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Fig. 7. - Component characteristics of trans- 
fer function for heptane vaporization 
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Fig. 8. - Comparison of heptane vaporization 
with burning rate characterized by n 
and i. 




